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_ • ' Abstract 
(N ' 

GRB 041006 was detected by HETE-2 at 12:18:08 UT on 06 October 2004. This GRB displays a soft 
X-ray emission, a precursor before the onset of the main event, and also a soft X-ray tail after the end 
of the main peak. The light curves in four different energy bands display different features; At higher 
energy bands several peaks are seen in the light curve, while at lower energy bands a single broader bump 
dominates. It is expected that these different features are the result of a mixture of several components 
each of which has different energetics and variability. To reveal the nature of each component, we analysed 
the time resolved spectra and they are successfully resolved into several components. We also found that 
these components can be classified into two distinct classes; One is a component which has an exponential 
decay of E p with a characteristic timescale shorter than ~ 30 sec, and its spectrum is well represented by 
a broken power law function, which is frequently observed in many prompt GRB emissions, so it should 
have an internal-shock origin. Another is a component whose E p is almost unchanged with characteristic 
timescale longer than ~ 60 sec, and shows a very soft emission and slower variability. The spectrum of the 
soft component is characterized by either a broken power law or a black body spectrum. This component 
might originate from a relatively wider and lower velocity jet or a photosphere of the fireball. By assuming 
that the soft component is a thermal emission, the radiation radius is initially 4.4 x 10 6 km, which is a 
typical radius of a blue supergiant, and its expansion velocity is 2.4 x 10 5 km s _1 in the source frame. 
Key words: gamma-rays:busts — X-rays: bursts — X-rays: individual (GRB041006) 
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1. Introduction 

On October 6, 2004 the High Energy Transient 
Explorer 2 (HETE-2) detected a gamma-ray burst (GRB) 
with soft X-ray emission before the onset of the main 



event. Such soft emission, a precursor, is predicted in 
some of theoretical models. The fireball undergoes a tran- 
sition from an optically thick phase to an optically thin 
phase, and thermal radiation (the fireball precursor) may 
occur during this transition (B. Paczynsky 1986; Daigne 
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& Mochkovitch 2002). A precursor (progenitor precursor) 
may also be emitted by the interaction of the jet with 
the progenitor star (Ramirez-Ruiz et al. 2002; Waxman 
& Meszaros 2003). The external shock by the first rela- 
tivistic shell can also produce the non-thermal precursor 
(Umcda et al. 2005). 

Soft precursors are occasionally detected in long GRBs. 
The first detection was made by the GINGA satellite 
(GRB 900126; Murakami et al. 1991). In more recent 
observations, the BeppoSAX (e.g. GRB 011121; Piro et 
al. 2005), HETE2 (e.g. GRB 030329; Vanderspek et al. 
2004) and Swift (e.g GRB 050820A; Cenko et al. 2006, 
GRB 060124; Romano et al. 2006, GRB 061121; Page et 
al. 2007) satellites have also detected precursors. Lazzati 
2005 studied bright long BATSE GRB light curves and 
found that in 20% of the cases there is evidence for soft 
emission before the main event. 

The precursor is usually detected as a single pulse that 
is well separated in time from the main event, typically 
several seconds to hundreds of seconds. The precursor of 
GRB 041006 is not well separated from the main event and 
is likely to be continuously active during the whole prompt 
GRB phase. Such a long lasting soft component was also 
observed in GRB 030329 (Vanderspek et al. 2004). Vetere 
et al. 2006 found that for some of the GRBs detected by 
the BeppoSAX, there is a slowly varying soft component 
underlying the highly variable main event. Borgonovo et 
al. 2007 analyzed the light curves obtained by BATSE, 
Konus, and BeppoSAX, and found that the width of the 
auto-correlation function shows a remarkable bimodal dis- 
tribution in the rest-frame of the source. This result sug- 
gests that there exists a slowly varying soft component in 
some GRBs. The relation between the underlying soft X- 
ray component, the X-ray precursor, and the main event 
is still open to question. 

In this paper, we present the results of multiple compo- 
nent analysis of the time resolved spectra of GRB 041006. 
Throughout this paper the peak energies are in the ob- 
server's frame, and quoted errors are at 90% C.L., unless 
specified otherwise. 

2. Observation 

GRB 041006 was detected with the HETE FREGATE 
(Atteia et al. 2003) and the WXM (Shirasaki et al. 2003) 
instruments at 12:18:08 UT on 06 October 2004 (Galassi 
et al. 2004). The WXM flight software localized the burst 
in real time, resulting in a GCN Notice 42 seconds after 
the burst trigger. The prompt error region was a cir- 
cle of 14 arcminute radius (90% confidence) centered at 
00 h 54 m 54 s , +01°18'37" (J2000). Ground analyses of the 
burst data allowed the error region to be refined to a cir- 
cle of 5.0 arcminute radius (90% confidence) centered at 
00 h 54 m 53 s , +01°12'04" (J2000). 

1.4 hours after the trigger, the optical afterglow was 
found by Da Costa et al. 2004, and the redshift was first re- 
ported by Fugazza et al. 2004 and later confirmed by Price 
et al. 2004 to be z = 0.716. Follow-up observations were 
made at various observation sites (e.g. Urata et al. 2007). 



VLA observations were made but no radio sources were 
detected (Soderberg et al. 2004). The X-ray afterglow was 
found by Butler et al. 2005, and it exhibited a power law 
decay with a slope of —1.0 ±0.1. The X-ray spectrum 
was characterized by an absorbed power law model with 
a photon index of T = 1.9±0.2 and n H = (1.1 ± 0.5) x 10 21 
cm~ 2 . The emergence of a supernova component was re- 
ported by Bikmaev et al. 2004 and Garg et al. 2004. The 
field of GRB 041006 was imaged by Soderberg et al. 2006 
using the WFC of the ACS on-board HST, and they found 
a SN 1998bw-like supernova dimmed by ^0.3 magnitudes. 

3. Analysis 

The data obtained by the WXM and FREGATE in- 
struments were reduced and calibrated in the standard 
manner. We used WXM TAG data and FREGATE PH 
data. 

3.1. Temporal Properties 

Figure 1 shows the light curves of GRB 041006 in four 
energy bands with 0.5 sec time resolution. T50 and Tgo 
are measured for each energy band, and they are shown 
in tabic 1. 

The burst can be divided into four major intervals ac- 
cording to spectral features, and each major interval is 
divided into a few sub-intervals for time-resolved spec- 
tral analysis. The time intervals for each sub-interval are 
shown in table 2. In interval 1 soft emission showing no 
prominent activity above 40 keV occurs, then harder emis- 
sions follow in intervals 2 and 3. In interval 4, the hard 
emission almost disappears and only gradually decaying 
soft emission is present. 

We call the emission seen in interval 1 an X-ray precur- 
sor. The precursor shows a structured light curve in the 
lowest energy band (2-10 keV), which indicates that two 
emissions are occurring successively . In interval 2, two 
peaks are seen in the higher energy bands (> 40 keV). 
The time history of the hardness ratio also clearly shows 
the corresponding peaks. In the lowest energy bands (< 
10 keV), structured emission is not clearly seen. In in- 
terval 3, two harder peaks arc seen in the highest energy 
band (80-400 keV), and this structure is less distinct in 
the lower energy bands. The emission in interval 4, which 
we call an X-ray tail, shows no prominent structure. 

From the dissimilarity of the light curves in the four 
energy bands, it is inferred that the total emission is com- 
posed of several independent emissions which have differ- 
ent characteristic energies. For an example, two compo- 
nents which contribute to the precursor, four components 
seen as a peak in the energy bands 40-80 keV and 80- 
400 keV, and one broad soft component which constitutes 
the major part of the light curve in the lowest energy 
band. To investigate this hypothesis, we performed time 
resolved spectral analysis based on a multiple-component 
spectrum model. 
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3.2. Average Spectral Properties 

The joint spectral analysis of WXM and FREGATE 
data was performed using XSPEC v. 11. 3.1 (Arnaud 1996). 
The time integrated spectrum of GRB 041006 is approx- 
imately described by a broken power law function (fig- 
ure 2); the low energy photon index is a = 1.28 ±0.02, 
the high energy index is j3 = 2.14 ±0.07, the break en- 
ergy is E p = 22.5 ± 1.7 keV and the flux at 1 keV is K = 
4.25 ±0.15 cm _2 s _1 keV -1 , where the quoted errors are 
one sigma. The x 2 is 111.19 for 79 dof, and Null hypothe- 
sis probability is 0.0099, so the fit is not very good. From 
this fitting result, we obtained S x = (5.24 ±0.08) x 10~ 6 
= (7.13±0.12) x 10 6 ergs cm 2 , where Sx 
and iS 7 denote fluences in the 2-30 keV and 30-400 keV 
energy ranges and the error is 1 sigma. As the ratio of 
fluences is log(S x /S 7 ) = —0.13, the GRB can be classified 
as an X-ray Rich GRB (Sakamoto et al. 2005). 

The isotropic energy is calculated from: 

E iso = —± / E$dE (1) 

where Dl is the luminosity distance, <E> is the differential 
photon spectrum, the range of energy intergation is from 
1 keV to 10000 keV in the source frame. We obtained 
E iso = 2.54^0 35 x 10 52 ergs. In figure 3, the peak energy 
in the source frame E p src is plotted against the isotropic 
energy Ei so (the point labeled "Total"). The relation for 
GRB 041006 obtained from the one component fit is com- 
pletely outside the Amati relation (Amati 2006). 

Looking at the residual plot in the top panel of figure 2, 
an additional soft component is apparently seen around 
6 keV and a systematic excess is also seen around 50- 
100 keV. Thus the total spectrum was fitted by a super- 
position of multiple basic functions. As basic functions, 
we considered a broken power law and a black-body. 

For the broken power law model, we used the following 
function to estimate the peak energy flux directly: 

A(E)=K/E%(E/E p )- a , E<E P (2) 

K/E 2 p {E/E p )-P, E>E p 

The parameters a and /?, which are the lower and higher 
energy photon indices, are restricted to the range of - 
2.0-2.0 and 2.5-5.0, respectively. The initial value of the 
break energy E p of the bknp basic function is determined 
from the local excess of the residual between the single 
bknp model and the observed data. The restriction to 
the break energy E p is applied so that the parameter con- 
verges around the initial value. 

The results of the spectral fit for three three-component 
models are shown in table 3. For comparison the result of 
the two-component model and a fit by the Band function 
(Band et al. 1993) and a broken power law function are 
also shown in the table. The fitting parameters for the 
models bbody*2+bknp and bknp*3 are given in table 4. 

Akaike's Information Criterion (AIC) is calculated for 
each model. AIC (Akaike 1974) is a very widely used 
criterion to evaluate the goodness of the statistical model 
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from both the goodness of fit and the complexity of the 
model. AIC is defined by the following equation: 

AIC = n\n(?- S J+2k, (3) 

where n is the number of data points, k is the number 
of free parameters to be estimated, and \ 2 1S the resid- 
ual sum of squares from the estimated model. The AIC 
includes a penalty that is an increasing function of the 
number of estimated parameters; overfitting is discour- 
aged, and thus this method enables one to find the best 
model for the data, with minimum of free parameters. 
The model with the lower value of AIC is the one to be 
preferred. 

The most preferable model is bbody*2+bknp. The 
model name is given by an algebraic expression of the 
name of a basic model. The second most preferable model 
is bknp*3. The AIC values for the two models are 6.87 
and 8.47 respectively. 

The lowest AIC does not necessarily select the true 
model, and the degree of the preference is estimated by 
the AIC difference. The relation between the degree of the 
preference and the AIC difference (Ax), however, depends 
on n and the models to be compared. So we evaluate 
the confidence limit of the AIC difference by carrying out 
a Monte Carlo simulation. The Monte Carlo simulation 
was performed by using the f akeit command of XSPEC, 
which generated 1000 PHA samples based on the spectral 
model to be tested. For each PHA sample, a spectral fit 
was performed for both the tested model and the model 
which gave the lowest AIC, and the AIC difference was 
calculated. 

The left panel of figure 4 shows a sim- 
ulated distribution of the AIC difference 

A bknp*3 = AIC bknp*3- AIC bbody*2+bknp- Tnc 

simulation was performed with the model spectrum 
bknp*3; the model parameters were obtained from the fit 
to the observed total spectrum. For each simulated PHA 
sample, model fit was performed for both the bknp*3 
model and bbody*2+bknp, which is the most preferred 
model. From this result the 90% confidence limit for 
^bknp*3 i s estimated as 4.7, below which 90% of samples 
are included. The observed AIC difference for the model 
bknp*3 is 2.64, so the model is acceptable at 90% C.L. 
In the case of the Band model (right hand panel of 
figure 4), for 98% of the samples the AIC is smaller than 
the most preferred model bbody*2+bknp. The observed 
AIC difference is 13.68, so the Band model is rejected 
at higher than 98% C.L. All the three three-component 
models are acceptable at 90% C.L. The two-component 
model is rejected at 90% C.L. 

As the time averaged spectrum of GRB 041006 is well 
represented by a superposition of the three components, 
we examined the E PtSTC -E lso relation for each one. The 
Ei so calculated for a model bknp*3 are summarized in 
table 7. The Ei so calculated for a model bbody*2+bknp 
is also shown in the table for the high energy component. 
The result are compared with the other GRBs in figure 3. 
The components with E p > 40 keV (C) and E p ~ 20 keV 
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(B) are well within the Amati relation, and the component 
E p ^6 keV (A) is out of the 90% distribution width of the 
Amati relation. The log(S x / S 7 ) for the three components 
are —0.3 for the component C, 0.78 for the component B, 
and 0.76 for the component A; thus they are classified as 
XRR, XRF and XRF, respectively. 

3.3. Time Resolved Spectral Properties 

Time resolved spectral analysis was performed for 12 
independent time intervals, and also for some intermedi- 
ate intervals which overlap part of one or two adjacent 
intervals to trace the spectral evolution more closely. We 
applied multi-component models in the spectral fit, where 
the model spectrum is represented as a superposition of an 
arbitrary number of basic functions. The basic functions 
considered here are black body (bbody), broken power law 
(bknp), and a single power law function (pi). The XSPEC 
built-in model is used for bbody and pi, for which the 
XSPEC model names arc bbodyrad and powerlaw respec- 
tively. For the broken power law model, we used equa- 
tion 2. 

The fitting results for various combinations of basic 
functions are summarized in table 5. The fitting parame- 
ters for the lowest AIC model are shown in table 6. The 
model spectra giving the lowest AIC at each interval are 
shown in figures 6 and 7. The expected number of compo- 
nents constituting the total spectrum is inferred from the 
number of local excesses in the residual plot for the bknp 
model, and also from the light curves in the four energy 
bands. As an example, the case of interval 2c is shown in 
figure 5. The spectrum is fitted with a single broken power 
law function, and E p is determined as ~20 keV. Looking 
at the residual plot shown in the bottom of the figure, 
local excesses around 6 keV and 60 keV are seen. So the 
spectrum of interval 2c is expected to be constituted from 
three components which have peak energies of 6, 20, and 
60 keV. In the case of interval 2b at least four components 
are expected from the light curves. One is the precursor 
component seen in interval 1, which is expected to be 
present in interval 2 if it is extrapolated smoothly. Two 
components corresponding to the two peaks seen in the 
40-80 keV energy band and one component correspond- 
ing to the broad soft emission in the lowest energy band 
are also expected to be present. So up to four components 
are examined for interval 2b. 

The model selection is carried out by examining the 
AIC difference, and the 90% confidence limit of the AIC 
difference is calculated by performing a Monte Carlo sim- 
ulation. By this statistical examination, single component 
models considered here are rejected for most of the inter- 
vals. The single component model is accepted only for 
intervals la, 4a, and 4b. For the other intervals, the sin- 
gle component model considered here is rejected at 90% 
C.L. and the multi-component models are preferred. 

For most of the intervals, the null hypothesis probabil- 
ity is larger than 0.1. For interval 2b, however, the null 
hypothesis probability is at most 0.003. This is proba- 
bly because unknown systematic errors are present in the 
data. 



4. Discussion 

The optical afterglow light curve in the R band can 
be fitted by a broken power-law model with a break time 
f fc = 0.16±0.04days (Staneket al. 2005). Taking^ as a jet 
break time and assuming a homogeneous density profile 
around the GRB, the jet opening angle 8 is estimated from 
the following equation (Sari et al. 1999, Nava et al. 2006): 

3/8 / x 1/8 

-) (^-) , (4) 



1 



E 



iso,52 



where no is the ambient particle density in cm~ 3 , rj 7 
the radiation efficiency, and £^0,52 = E iso / (10 52 erg). 
Assuming hq = 3 and ?7 7 = 0.2, we obtain a jet open- 
ing angle of 3.4°. If the GRB is viewed on-axis, the 
collimation-corrected total energy can be estimated from 
E 7 = (1 — cos 6) Ei SO . The corrected total energies for 
the three components are 2.4±5'.I° x 10 49 erg for E^src = 



123If? keV, 0.49Io.2 x 10 erg for E p ^ c = 44+^ keV, 
and 1.7±f* x 10 48 erg for £ p , src = 8.4±?;g keV. These val- 
ues do not follow the Ghirlanda relation (Ghirlanda et al. 
2007) except for the component with E p ~ 6 keV. That is, 
the -Ep jS rc expected from the Ghirlanda relation are 39.4, 
13.0 and 6.2 keV for the components with E p > 40 keV, 
~ 20 keV, and ~ 6 keV, respectively Taking a 5% un- 
certainty in the Ghirlanda relation, the observed E p for 
the the components with E p > 40 keV and ~ 20 keV are 
incompatible. 

We also tested the Liang-Zhang relation (Liang & 
Zhang 2005). The isotropic energies E- lSQ ^2 calculated 
by equation (5) of Liang & Zhang 2005 are: 2.54, 0.132, 
3.28 and 24.1 for components "total", A, B, and C, re- 
spectively. The isotropic energy derived from the fit to 
a single broken power law function are consistent with 
the isotropic energy derived from the Liang- Zhang rela- 
tion. On the other hand, the isotropic energies derived 
for components B and C are incompatible with those ob- 
tained from the relation. 

Looking at the time evolution of E p obtained by the 
time resolved spectral analysis shown in figure 8, we can 
identify seven components. Each component is interpo- 
lated with a solid line, and is given an identifier A, Bi, 
B 2 , Ci, C 2 , C 3 or C 4 . 

The most preferred spectral model for component A in 
interval la is the bbody model. The calculated emission 
radius is 4.351*' 4 x 10 6 km, which corresponds to 6 solar 
radii and is a typical radius for a blue supergiant. The 
AIC difference for the second-most preferred bknp model 
is 3.31 and its 90% confidence limit is 4.9, so the bknp is 
also acceptable. The AIC differences for the power law 
spectrum with and without absorption (wabs*pl and pi) 
are larger than 8.9, and their 90% confidence limits are 
less than 0.3, so these models are rejected at 90% C.L. 

For interval lb, the acceptable models arc bbody*2, 
bbody+bknp and bknp*2, all of which are two-component 
models. None of the single component models considered 
here is preferable and all are rejected at 90% C.L. Thus 
it is likely that the emission in interval lb is composed of 
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two components (A and Bi). The spectral type of each 
component is not uniquely determined from this result; 
it is either a black body or a broken power law function. 
Assuming that component Bi is black body radiation, the 
calculated emission radius is about one solar radius. 

In intervals 2a-2d, the soft components A and Bi are 
present in all the acceptable models. The peak energies of 
the components are almost constant during intervals 1 and 
2, and they decrease slowly, with decay time 72±42 sec for 
component A and 57±33 sec for component Bi. Assuming 
that the components originate from thermal emission, we 
can derive the evolution of the radiation radii, and they 
are shown in figure 9 with the filled circles for compo- 
nent A and with open circles for component B±. The 
data points for component B\ are shifted by a factor of 
four. The data points for intervals 1 and 2 are fitted with 
a linear function, and we calculate the apparent expan- 
sion velocity for component A to be (6.3 ± 1.5) x 10 5 km 
s _1 , which is twice the speed of light. This supcrlumi- 
nal motion is observed when the emitter is moving with 
relativistic velocity toward the observer. The relation be- 
tween the apparent expansion velocity v and the velocity 
measured in the source frame v' is given by: 

v= -. (5) 

(l + *)(l-£) 

The expansion velocity in the source frame is 2.35 x 10 5 km 
s _1 , and the corresponding Lorenz factor is 1.6. The ap- 
parent expansion rate for component Bi is found to be 

1.1 x 10 5 km s _1 , and the velocity in the source frame is 

1.2 x 10 5 km s _1 , which is half the velocity of component 
A. This result indicates that the soft component originates 
from the GRB photosphere expanding with a mildly rel- 
ativistic speed. According to the current models of GRB 
photosphere (e.g. Meszaros et al. 2002; Rces & Meszaros 
2005), however, it is difficult to interpret a blackbody with 
essentially the same temperature but an increasing radius, 
unless the temperature is boosted by the growing Lorentz 
factor of the photosphere. 

If the component originates in an internal shock accord- 
ing to the model of Zhang & Meszaros 2002 the following 
relation should be satisfied: 

E p (x L 1/2 T- 2 (6) 

where L is the luminosity and T is the bulk Lorentz factor 
of the shock. If the spectral shape does not change, the 
normalization constant K of equation 2 is proportional to 
the luminosity. As the a and (3 are not well constrained 
in the multi-component model due to the correlation of 
the parameters among the components, the luminosity is 
not well constrained. We have plotted the E p -K relation 
in figure 10. If V is constant and the spectral shape does 
not change during the emission, we expect that E p will 
be proportional to K x l 2 . No clear correlation is found 
for component A (filled circle). For component B\ (filled 
triangle) the expected correlation is not found either, and 
it shows a negative correlation. 

The higher energy components of the interval 2, Ci and 
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C2, which correspond to the two peaks seen in the 40- 
80 kcV light curve, are resolved as a broken power law 
spectrum for which E p is around 50-90 keV. If we assume 
that E p decreases exponentially as seen in many GRBs, we 
can derive the correspondence among the E p as indicated 
in figure 8. The decay constant of the E p is ~20 sec. 

At interval 3, the first precursor component seen in in- 
terval la (component A) is not well resolved. Component 
B2 has a similar E p to that of component Bi, but its E p is 
somehow systematically higher than the extrapolation of 
Bi. Assuming that B2 is thermal emission, its radiation 
radius is calculated and shown in figure 9. The radiation 
radius is well below the extrapolation of those for Bi. The 
Ep-K relation of B2 is shown in figure 10, and it does not 
follow the relation given by equation 6. 

The highly variable spectra whose emission peaks vary 
from 100 keV to 40 keV are also resolved (C3, C4), and 
they correspond to the emissions seen in the light curve 
of the highest energy band. From figure 8, the E p of the 
components decrease exponentially with time with a decay 
constant of ^5 sec. 

The Ep-K relations for components Ci, C2, C3 and C4 
are also shown in figure 10. Although there are few data 
points for each component, the E p -K relation is satisfied 
except for two points. Both the exceptions are at the time 
intervals corresponding to the rising part of the compo- 
nents Ci and C3. During the rise, due to the curvature 
effect, the emission from a part of the shock front that 
is moving toward us dominates. After that, the emission 
is averaged over a wider region, so the emission proper- 
ties may change between the rising part and the following 
part. 

In interval 4a, component B2 is likely to remain and a 
black body spectrum with T = 1 keV or a broken power 
law spectrum with E p ~ 4 keV is also likely to be present. 
In interval 4b, a power law spectrum with photon index 
1.9 is the most preferred model, which is almost the same 
as the afterglow spectrum observed by Chandra. 

5. Conclusion 

We have analyzed the time resolved spectra of 
GRB 041006 and successfully resolved the components 
corresponding to the hard spikes and the soft broad bump 
observed in the multi-energy band light curves. The com- 
ponents may be divided into two classes. One is compo- 
nent A, which has almost constant E p around 6 keV, and 
components Bi and B2 which have almost constant E p 
around 20 keV. E p for this class gradually decreases on a 
timescale, 60-70 s. The spectral type is well represented 
by a broken power law function or a black body radiation 
function. Assuming that the emission of this component is 
due to black body radiation, we derived the emission radii. 
At the beginning of the emission they are 4xl0 6 km for 
component A and 7xl0 5 km for components Bi and B2. 
The expansion velocity in the source frame is also derived; 
it is 0.78 c and 0.4 c for components A and Bi, respec- 
tively. The emission radius of component B2 is almost 
constant. 
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Table 1. Temporal properties, T 50 and Tg , of GRB 041006. 



energy range 


T50* 


T90* 




(kcV) 


00 


(s) 


2 


- 10 


13.9±0.08 


38.2±0.40 


10 


- 25 


11.9±0.16 


27.3±1.44 


40 


- 80 


10.2±0.09 


19.6±0.10 


80 


- 400 


3.7 ±0.25 


17.4±0.25 



* The quoted errors correspond to one sigma. 



The i^-Luminosity relation is examined for these com- 
ponents and compared with the prediction of the internal 
shock model. We used a normalization constant K in 
equation 2 instead of deriving the luminosity. According 
to the internal shock model of Zhang & Meszaros 2002. 
E p is proportional to L 1 / 2 if the bulk Lorentz factor of the 
shock is constant during the emission. We could not find 
such a correlation for components A, Bi and B2. 

The second class comprises the components whose E p is 
larger than the former class and shows a relatively rapid 
decrease on a timescale of 5-20 sec. The spectra are well 
represented by a broken power law function, and the E p -K 
relation almost follows the relation expected for an inter- 
nal shock origin, so this could explain their origin. 

We could not reach a conclusion about the origin of 
the soft component observed for GRB 041006. However, 
the difference in its time variability with respect to the 
higher energy component suggests that it originates from 
different emission sites, such as acceleration by a wider jet, 
emission from a supernova shock breakout, or emission 
from the photosphere of the fireball. 



Table 2. Time intervals used for 

time resolved spectral analysis. 

time interval start* - end* 



(s) 



la 


2.5 - 6.0 


11) 


6.0 - 12.5 


2a 


12.5 - 16.5 


21) 


16.5 - 19.5 


2c 


19.5 - 23.0 


2d 


23.0 - 27.5 


3 a 


27.5 - 29.5 


31) 


29.5 - 31.0 


3c 


31.0 - 34.0 


3d 


34.0 - 38.0 


4a 


38.0 - 42.5 


4b 


42.5 - 60.0 


2a' 


15.0 - 16.5 


2c" 


22.0 - 24.0 


3b' 


30.0 - 32.0 


3c" 


33.0 - 35.0 



* The offset time is the trigger time 2004-10-06 12:18:08.63933. 
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Table 7. Isotropic energies Ej so 52 and rest-frame peak 



energies E p , Brc derived from the 


average spectrum. 


Component 


Ep^src 


-Biso,52 




(keV) 




Total* 


38.6 ± 2.9 


/ - 04 -0.35 


At 




0.0941°$ 


Bt 




0.28+0? 


ct 


123 +2 ? 


1 36+ 04 

i " 3D -0.8 


C't 


126+g 


1 09+O.5 
J - oz — 0.3 


* bknp model, 
t bknp*3 model. 



bbody*2+bknp model. 
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Table 3. Results of the spectral fit to the time averaged spectrum. 



model 


n* 


fct 


x n 






A x # 


parameters** 


bbody*2+bknp 


83 


8 


74.35 


0.499 


6.87 




T= 


■■1A,5.5,E P =74 


bknp*3 


83 


12 


68.84 


0.551 


8.47 


1.6(4.7) 


E p - 


=5,25,72 


bbody+bknp*2 
bknp* 2 


83 


10 


73.75 


0.453 


10.19 


3.32(4.1) 


T= 


:1.6,£!p=23,73 


83 


8 


77.80 


0.390 


10.63 


3.76(<0) 


E P 


=5,24 


band 


83 


4 


96.55 


0.087 


20.55 


13.68(<0) 


E p 


=38 


bknp 


83 


4 


111.19 


0.010 


32.27 


25.40(<0) 


E p 


=22 



" IN umbers ol data points used tor the lit. 

t Numbers of model parameters. 

i Chi-square of the fit 

§ Null hypothesis probability. 

II Akaikc information criterion. 

# AIC difference between the corresponding model and the lowest AIC model. The numbers in parentheses represent the 
90% confidence limits of the AIC difference. 

** T is the black body temperature in keV and E p is the break energy of the brknp model in keV. 



Table 4. Fitting parameters for the time averaged spectrum. 



model 


component 


parameters* 












bbody*2+bknp 


1 


T= 1.40 +, (ji6 


-^bbody = 0.16 ±0.04 












2 
3 


nip — *o.o_ 15 6 


A'bbo dy = 0.44 ±0.10 
a = 1.33±°?4 


= 






= 37.8±|:? 


bknp*3 


1 

2 
3 


E p = 71.9_n g 

Ep = 25.4±|;g 

F — 4 Q+ 13 
£jp — *.J_ 06 


a = 1.3^| 

a = 1.2153 
a = -2.00± 3 ;° 


fS = 
P = 
(3 = 


s.oot™ 


-^bknp 
-Kbknp 
^ bknp 


= 43.4+^ 
= 19-8^ 4 3 
= 3- 69 ti:o 



* T and ftTbbody = ^km/^io arc * nc temperature in units of keV and normalization constant for the black-body radiation 
model, respectively. /?km is the source radius in units of km. Dio is the distance to the source in units of 10 kpc. Ep, 
a, (3, -Kbknp are the break energy in units of keV, low energy photon index, high energy photon index, and normalization 
constant defined in equation 2. -R"bknp is in units of keV cm" 2 s — 1 . 
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Table 5. Results of spectral model fitting to the time resolved spectra. 



in fpwfi 1 


rn < » r ol 

1.L1U l_X<^± 




fct 


Y 2t 
A 


7^§ 
P 


Arc 11 




l^CU. (TX 41j1Lj Lt^jJ. O 


1 a 


r» ri fin^r 


52 

■> z 


2 


41 38 


809 
u .ouz 


-7 87 

i .O I 




T=2 




bknp 


52 

UZi 


4 


40 75 


769 

U . f uz 


-4 68 


O . i- c/\o. Cf J 


R —7 3 

-L/p 1 .O 




wabs*pl 


52 


3 


47.26 


0.544 


1.03 


8.90(1.1) 


Ct~3.0.7lH = 16 




d1 


52 


2 


56.57 


0.243 


8.38 


16 25f<0) 


a=2.1 


1 h 

X 1 J 


uuut-iy z 


^9 


4 


97 


SQ^ 

u.oyo 


1 73 

-LU. 1 O 




7^_1 /ICQ 






52 


(j 


35.92 


0.857 


-7.24 


3.49(4.2) 


J. -L . O j i-J'p OU 




bknr>*9 

UlYilU Z 


5 r > 

u z 


s 

o 


35 60 

oo .UU 


813 


-3 70 

. 1 U 


7 03f7 4~) 

f . UO 1 i . L ±j 


-1— /p — u,ou 




bknp 


52 


4 


42.92 


0.681 


-1.98 


8.75(<0) 


-L^P U 




bboH v-4-r>l 


52 


4 


49.93 


0.396 


5.89 


16 62r<0) 


T=2.1,a=1.9 




d1 


52 


2 


63.52 


0.095 


14.41 


25.14(<0) 


p=1.9 


9a 

Ztl 


bbndv*9+hknn 


SO 

O V.' 


a 

O 


59 34 


857 

U .OO i 
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f . C/U 




T— 1 7 5 9 F, —84 
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UxvilU Z 


80 

o u 
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O 


61.24 


813 
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o .OO 
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U U KJ KA. y [ L J'XYJ.i L/ Zj 
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10 
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0.837 
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7^_o „ r — 2S 83 

J. Z* . U i J-Jp Zj o , o o 




UxVllJJ O 


Sfl 


1 9 
r ^ 


^7 fiS 


810 

U.OIU 


9 1 7 


^ 73(^0 4~> 
o . i o i y .4 t 


X^p O ,Zi4:,00 




bknp 


80 


4 


70.48 


0.657 


-2.13 


5.77(0.5) 


i? p =25 


2b 


bbody*2+bknp 


80 


g 


104.91 


0.007 


37.69 




T— 1 4 5 4 K —84 
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12 


99 33 


008 
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1 01 08 
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4:. 1 O I O . Z I 
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bbody+bknp*2 


80 


10 


106.05 


0.004 


42.55 


5.22(5.5) 


r=1.5,£?p=22,85 


9c 


bbndv*9+hknn*9 

U U KJ Kl y Z \^ UxViX VJ Z 

bbodv*2-4-bknn 

UUU KXy Z UxYll 


73 


12 


49 53 


853 

u .ooo 


-4.32 




T— 1 3 5 B —59 98 

J- — J. . O . O . U j J-/p — O Zj , o 




73 


8 


56.66 


0.760 


-2.50 


1.67(<0) 


T— 1 3 5 F —53 

-L - — -L . O j O . U j J-Jp O O 




V» l~i n rl t 7 4- i~» Irn r\ ^ 9 
UUU IX y U ±\.iX|J Z 


1 '_> 


i ii 

1 u 


^fi fii 

UU.Ul 


709 

U. lUi 


1 44 


^ 7fi^0 9^ 
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UlYilU O 


73 


12 


53 58 


739 

U . f OtJ 


1.42 


5 74f0 9^ 

O . 1 t:1 U. Z J 


F — 5 5 1 8 74 




bknp* 2 


73 


8 


62.24 


0.574 


4.36 


8 68f0 06^ 

O.UUI U.UU / 


B — 19 54 




bbody-f-bknp 


73 


G 


66.70 


0.488 


5.41 


9.73(<0) 


T=4.7.£'p=55 




bknp 


73 


4 


87.99 


0.006 


21.63 


25.72(<0) 


Bp=23 


9d 


hhnHv*9-l-hVnn 

UUUUV Z ulYil jJ 


fifi 

uu 


a 
o 


fi/L 70 

U^. ( u 


9^.4 


1 4 RQ 

±4:. Uty 




T— 1 9 4 fi B — R9 
j — r . Zi ; 4.u — uZi 




rt ri f 1 n v -4- n k lTn 

UUU *-lY 1^ L *J\.I.l JJ 


66 


(j 


72.12 


0.136 


17.85 


3.16(0.9) 


T— 4 5 B —62 

J_ t: . U j J-Jp U Zj 




*hlnrn*9 

UxVllJJ Z 


fifi 

uu 


a 
o 


7D ^9 
( u .oo 


1 9Q 


90 1 Q 
^u. i y 


^ ^on 9^ 


p -inn 
j-jp — io,<jy 




bknp 


GG 


4 


80.21 


0.060 


20.87 


6.18(0.1) 


Sp=18 




bbody+bknp*2 


GG 


10 


67.50 


0.140 


21.48 


6.79(5.5) 


T— 1 6 B —17 60 




hknn*3 

UivllU rj 


GG 


12 


66.84 


0.113 


24.83 


10.14(4.4) 


B —4 1 7 60 

J-Jp 1 I .UU 
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53 
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74 


o 


63 79 
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U.OO I 


4 93 


4.40(4.9) 
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74 


10 


61.83 
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6.18(6.8) 


T—p. n B —50 92 
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bknp 


74 


4 


75.48 
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74 
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1 f,9(~\ 1 8"! 
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8 
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1 u 


fiQ 1 Q 


AXl 

U .4:0 / 


1 ^ fiQ 


3 ^8^8 fi^ 

O.UOl O.U 1 


B — fi 10 91 84 
J-/p — u, r u , ^ r ,04 




hhnHv+hVrm*9 

U UUlXy T UxyIIlU 


84 


1 


80 1 7 
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3 97f4 0^ 
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J. — O^J-Jp — OO 




bknp*3 


84 


12 


79.88 


0.245 


19.78 


7.67(7.2) 


£p=5,26,80 




bknp 


84 


4 


107.35 


0.022 


28.60 


16.49(<0) 


E P =Q7 


3c 


bknp*3 


73 


12 


70.36 


0.193 


21.32 




S p =26,44,120 




bbody+bknp*3 


73 


14 


67.43 


0.211 


22.20 


0.88(4.5) 


7=1.2,2^=26,44,118 




bknp*2 


73 


8 


80.75 


0.090 


23.37 


2.05(2.3) 


£; p =44,130 




bbody+bknp*2 


73 


10 


78.07 


0.096 


24.90 


3.58(1.3) 


T=l.l,2? p =44,117 




bknp*4 


73 


16 


67.91 


0.153 


26.72 


5.40(7.4) 


r5 p =6,26,44,119 




bknp 


73 


4 


98.92 


0.011 


30.18 


8.86(<0) 


r5 p =56 


3d 


bbody+bknp 


80 


6 


76.28 


0.405 


8.19 




7=6.1,2^=72 
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interval 


model 


n* 


fct 


x n 


p§ 


AIC" 


A x # 


parameters** 




bknp* 2 


80 


8 


77.40 


0.310 


13.36 


5.17(5.8) 


Ep=21,47 




bknp 


80 


4 


86.42 


0.194 


14.18 


5.99(<0) 


£ p =24 




bknp*3 


80 


12 


74.91 


0.264 


18.74 


10.55(13.6) 


£ p =23,43,75 


4a 


bbody*2 


66 


4 


59.23 


0.576 


0.86 




T=1.2,5.2 




bbody+bknp 


66 


6 


59.14 


0.505 


4.76 


3.90(7.1) 


T=1.2,E p =2A 




bknp 


66 


4 


63.09 


0.438 


5.02 


4.16(2.8) 


E p =26 




bknp* 2 


66 


8 


57.36 


0.496 


6.74 


5.88(7.4) 


E p =4,25 




bbody+pl 


GG 


4 


73.06 


0.159 


14.71 


13.85(1.4) 


T=4.7,a=2.3 




Pi 


66 


2 


100.05 


0.003 


31.46 


30.60(<0) 


a=2.0 



4b 



pl 

bbody+pl 

bknp 

bbody 



52 2 

52 4 

52 4 

52 2 



47.31 0.582 

44.82 0.604 

45.13 0.591 

69.71 0.034 



-0.92 
0.27 
0.63 

19.24 



1.19(3.1) 
1.55(3.6) 
20.16(<0) 



a=1.9 
T=1.5,a= 



T= 



=4 
1.7 



" JN umbers ot data points used tor the tit. 

* Numbers of model parameters, 
i Chi-square of the fit 

§ Null hypothesis probability. 
II Akaike information criterion. 

# AIC difference between the corresponding model and the lowest AIC model. The numbers in parentheses represent the 
90% confidence limits of the AIC difference. 

** T is the black body temperature in units of keV, E p is the break energy of the brknp model in units of keV, a is the 
power law photon index of the pl model, and njj is the column density measured in unit 10 22 . 
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Table 6. Fitting parameters for the most preferred models, that is, the model that gives the lowest AIC. 



lTlt"PT"\r£1 1 
llluvl Vdl 


pAmnnn r*n1~ 

\J _L± 1 IJ KJ 1 1L _L± Li 
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T — 1 Q9+ uau 
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-* v bbody — 4: *- L ' —1.4 A lu 
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-Kbbody = 3.95^33 


= 
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+5.8 
-11 
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I 
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3 


t — i 4n+ u ' 2a 

T = 5.40l°;g 
E — 84 3+™ 4 
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-"-bbody — 1 - uz — 0.63 u 
-Kbbody = 13.0l 5 ' 7 
a - 1 26 + S1o 


— 


5 00+°,-°.°, 


Kui — 57 8 

1 bknp u ' -° 


+13.9 
-12.2 


2c 


1 

2 
3 
4 


J — i>ofl -0.077 

T — E5 01+iL 
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p — T . 6 

= 95.5 g y 


/Vbbody = 1.44™Xl0 3 

Kuu a - 25 + ?o 9 
a = 0.24+1;° 
a = O.OGl^ 4 


= 

0= 
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5.ooi?;°° 


Kuv^r* = 97 9 

' 1 DKnp " 1 

A'bknp = 78.4 


+3 5 
-40 
+ 19 
-50 
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1 

2 
3 


r = i.2C£S 

T = 4.65±g;g 
^ = 62.1^ 


^bbody = 1.01^Xl0 3 
-Kbbody = 26.3l 9 ' 4 

a = 1.22+°;? 


0= 


5.00l°;° 


-Kbknp = 54.1 


+ 11.5 
-10.9 


3a 


1 

2 


T = 6.8tJ'j 
E p = 95.8+^ 5 5 


-^bbody = 3.6ll 2 '5 
1 trn+0.07 

a = l- 50 -o.o7 


0= 


5.00+°-° 


i^bknp = 107^ 


-IT 

-18 


3b 


1 

2 


-Ep — 25.3 2 q 

E p = 81.9l 9 ; 9 


a= -0.92+J'J 
« = 1.05±°;^ 


0= 

= 


5.00l°;° 

q oq+0.52 
°' ZO -0.46 


K = 68.7±£ 
K = 386l^ 2 




3c 


1 

2 
3 


E p = 25.8+4Q 
Ep = 44.0i 3 3 6 
Ep = 119±£ 


a = -0.10+°f 

a = - 2 -00±o.oo 
— 1.33 q 


= 

= 
= 


5.00+°° 

2-66l^;° 9 

5.ool?;°° 


K = 68.1+^ 
71=115+3° 
K = 159l°g 




3d 


1 

2 


T = 6.051^ 
E -71 9+ 14 


-fr'bbody = 5. 18li g 

_ i 39+0.05 


= 


4.32l°J 8 


K = 55.7l r 2 




4a 


1 

2 


7 1 - 1 9Q+»-i« 
— 1 - zo -0.16 

T = 5.161°;^ 


/Vbbody = 8.09±§-; x 10 2 

^bbody = 4.66±f"g 










4b 


1 


OL — 1.93_q _-[_4 


A pl - 2 - 74 -0.68 











* T and A" bbo( j y = R^ m /D'± are the temperature in units of keV and normalization constant for the black-body radiation 
model, respectively, ii^m is the source radius in units of km. Dio is the distance to the source in units of 10 kpc. E p , 



a, p, Kbknp are the break energy in units of keV, low energy photon index, high energy photon index, and normalization 
constant defined in equation 2. The unit of i^bknp i s keV cm~ 2 s" 1 . K p \ is the normalization constant for power law 
spectrum defined as photon flux at 1 keV in unit of photons keV -1 cm -2 s _1 . 
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Fig. 1. Light curves of GRB041006 in four energy bands and hardness ratio. From top to bottom, 2—10 keV, 10-25 keV, 
40-80 keV, and 80-400 keV. The hardness ratio is calculated by dividing the 40-80 keV count rate by the 2-10 keV 
count rate. The vertical lines represent the boundaries of the time intervals for time resolved spectral analysis. 
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Fig. 2. Time averaged unfolded spectrum expressed in vf u . Top: Fitting result for the broken power law model. Bottom: Fitting 
result for the three-component model represented by a superposition of one broken power law function and two blackbody functions. 
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Fig. 3. i?p,src--Eiso relation for long GRBs. The open circles represent the long GRBs compiled by Amati 2006. The 
solid circles represent GRB 041006. The solid circle labeled "Total" is derived from a single broken power law model. 
The other solid circles are derived from a bknp*3 model. The parameters obtained in this work are summarized in ta- 
ble 7. The solid line represent the average relation derived from all the points of the open circles, while the dashed 
lines represent lower and upper boundaries which arc parallel to the average relation and contain 90% of the points. 




AIC(bknp*3) - AIC(bbody*2+bknp) AIC(Band) - AIC(bbody*2+bknp) 



Fig. 4. Left: Simulated distribution of AIC differences between the bknp*3 and bbody*2+bknp models (left). The simulation is 
performed using the bknp*3 model, and model fitting to the simulated data is carried out for both the models. The bbody*2+bknp 
model is the most preferable model for the time integrated spectrum. The AICs for the two models are calculated for each simulated 
spectrum. The fraction of events with Aaic > corresponds to the probability of selecting the wrong model. Right: Same plot for 
the Band model. 
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Fig. 5. An example of spectral fitting for interval 2c, where a single-component model is used. 
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Fig. 6. Time resolved unfolded spectra for intervals 1 and 2. The residual between the observation and the model is also shown at 
the bottom panel of each figure. The spectrum is expressed in vj v . The most preferable model spectra are plotted as a solid line 
(total) and dashed lines (basic function). 
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Fig. 7. Time resolved unfolded spectra for interval 3 and 4. 
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Fig. 8. Peak energy calculated for each interval by fitting the data with multi-component models. The points which are 
inferred to belong to identical components are interpolated with a line. The vertical error bar corresponds to 90% C.L. 



18 



Y. Shirasaki ct al. 



[Vol. , 



25 



20 



- 1 1 r~ 



E 



C 15 
w 

T3 



o 

H — ' 

CO 

T3 
CO 

cn 



10 



• Component A 

o Component B1 and B2 (x 4) 




10 20 30 

Time since the trigger (s) 



40 



Fig. 9. Evolution of the radiation radii of the black body components. The rilled circles represent compo- 
nent A of figure 8. The open circles represent components Bi and B2, for which the radius is mul- 
tiplied by four. The solid and dashed lines represent the linear fit to the data of intervals 1 and 2. 




Fig. 10. The relation between E p and K of equation 2 for each component. Solid lines represent the relation E p oc K' 
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